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There are three regimes of gravitational-radiation-reaction-induced inspiral for a compact body with 
mass u, in a circular, equatorial orbit around a Kerr black hole with mass M ^> /i: (i) The adiabatic 
inspiral regime, in which the body gradually descends through a sequence of circular, geodesic orbits, 
(ii) A transition regime, near the innermost stable circular orbit (isco). (iii) The plunge regime, in 
which the body travels on a geodesic from slightly below the isco into the hole's horizon. This paper 
gives an analytic treatment of the transition regime and shows that, with some luck, gravitational 
waves from the transition might be measurable by the space-based LISA mission. 
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PACS numbers: 04.30.Db, 04.80. Nn, 97.60.Lf 



I. INTRODUCTION AND SUMMARY 

The space-based Laser Interferometer Space Antenna 
(LISA) Q, if it flies, is likely to detect and study the 
gravitational waves from white dwarfs, neutron stars and 
small black holes with masses fj, > 1M©, spiraling into 
Massive (M ~ 1O 5 -1O 8 M > fi) black holes in the nuclei 
of distant galaxies |^ [|]. In preparation for these studies, 
it is necessary to understand, theoretically, the radiation- 
reaction-induced evolution of the inspiral orbits, and the 
gravitational waveforms that they emit. 

Regardless of an orbit's shape and orientation, when 
H <C M the orbital evolution can be divided into three 
regimes: (i) The adiabatic inspiral regime, in which the 
body gradually descends through a sequence of geodesic 
orbits with gradually changing "constants" of the mo- 
tion E = (energy), L — (polar component of angular 
momentum), and Q — (Carter constant), (ii) A transi- 
tion regime, in which the character of the orbit gradually 
changes from inspiral to plunge, (iii) A plunge regime, in 
which the body plunges into the horizon along a geodesic 
with (nearly) unchanging E, L and Q. 

The plunge regime, being (essentially) ordinary 
geodesic motion, is well understood; and the adiabatic 
inspiral regime is the focus of extensive current research 
(see, e.g. [0-p[). By contrast, so far as we are aware, 
there have been no publications dealing with the transi- 
tion regime. 

In 1990-91, we carried out an initial exploration of the 
transition regime for the special case of circular, equa- 
torial orbits; but due to the press of other projects we 
did not publish it. Now, with prospects for LISA looking 
good, we have resurrected our work and present it here, 
in the context of LISA. 

We begin, in Sec. Q, by summarizing some key, well- 
known detai ls of the inspiral and plunge regimes. Then 
in Sec. Ill A we present a qualitative picture of the tran- 
sition from inspiral to plunge, based on the motion of a 



particle in a slowly changing effective potential (Fig. jl|). 
With the aid of this qualitative picture, in Sec. Ill B| we 
derive a non-geodesic equ ation of motion for the tran- 
sition regime, and in Sec. Ill C we construct the solu- 
tion to tha t equation of motion (Figs. |^ and ||). Then 



in Sec. [V, with the aid of our solution, we estimate 
the gravitational- wave signal strength from the transition 
regime and the signal-to-noise ratio that it would produce 
in LISA. We conclude that, with some luck, LISA may be 
able to detect and study the transition waves. In Sec. 
we make concluding remarks about the need for further 
research. 



II. ADIABATIC INSPIRAL AND PLUNGE 



Throughout this paper we use Boyer-Lindquist coor- 
dinates (t, r, 0, 4>) iQ for the massive hole's Kerr met- 
ric, and we use geometrized units, with G = c = 1. 
The hole's mass is M and the inspiraling body's mass 
is /i = r\M . We use M and /i to construct dimension- 
less versions (denoted by tildes) of many dimensionfull 
quantities; for example, f — r/M, and t — t/M. The 
hole's dimcnsionlcss spin parameter is a = (spin angular 
momentum) /M 2 (with —1 < a < +1). The body moves 
around its circular, equatorial orbit in the +0 direction, 
so a > corresponds to an orbit that is prograde relative 
to the hole's spin, and a < to a retrograde orbit. 

When the inspiraling body is not too close to the in- 
nermost stable circular orbit (isco), it moves on a circular 
geodesic orbit with dimensionless angular velocity M 



alt r 3 / 2 
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(where cj> is angle around the orbit) and with orbital en- 
ergy u 

1 - 2/f + a/f 3 / 2 
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As it moves, the body radiates energy into gravita- 
tional waves at a rate given by M 



5 



(2.3) 



where £ is a general relativistic correction to the Newto- 
nian, quadrupole- moment formula (Table II of Ref. Q). 
This energy loss causes the orbit to shrink adiabatically 
at a rate given by 



dr —E, 



GW 



dt dE/dr 



(2.4) 



The inspiral continues adiabatically until the body 
nears the isco, which is at the dimensionless radius 
risco = r isco /M given by @] 

f isco = 3 + Z 2 - sign(o)[(3 - Z x )(3 + Z x + 2Z 2 )] X / 2 , 

Zi = 1 + (1 - a 2 ) 1/3 [(l + «) V3 + (1 - «) V3 ] , 
Z 2 = (3a 2 + Zx 2 ) 1 / 2 ; (2.5) 

cf. Table |. The circular geodesic orbit at the isco has 
dimensionless angular velocity (Table |]), energy, and an- 
gular momentum given by j^J^] 



-3/2 . 
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As the body nears the isco, its inspiral gradually ceases 
to be adiabatic and it enters the transition regime (Sec. 
III). Radiation reaction (as controlled by Eqw) contin- 
ues to drive the orbital evolution throughout the transi- 
tion regime, but gradually becomes unimportant as the 
transition ends and pure plunge takes over. 

The plunge is described to high accuracy by reaction- 
free geodesic motion; Eqs. (33.32) of Ref. 0. Up to frac- 
tional corrections of order rf/ b , the orbital energy and 
angular momentum of the plunging body are eq ual to 
-Eisco and £i sco throughout the plunge, [cf. Eq. (3.2£) 
below] . 



III. THE TRANSITION FROM ADIABATIC 
INSPIRAL TO PLUNGE 

A. Qualitative Explanation of Transition 

As the body nears its innermost stable circular orbit, 
T = ^isco, the adiabatic approximation begins to break 



down. This breakdown can be understood in terms of the 
effective potential, which governs geodesic radial motion 
via the equation 



/df 



~\ 2 



dr\ - ~ - 

|:J =E 2 ~V(?,E,L) , (3.1) 



where E = E/fx = E/(r]M), L = L/(fj,M) = L/(r]M 2 ), 
and f = t/M are the body's dimensionless energy, an- 
gular momentum, and proper time. The explicit form of 
the effective potential can be inferred from Eqs. (33.32) 
and (33.33) of MTW 0: 



V{f,E,L)=E 2 - 



+ a 2 ) - La] 



2f + a 2 )[r 2 + (L - Ea) 2 }) . (3.2) 
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For a Schwarzschild black hole, this reduces to 



V(f,E,L)= (l 




for a = (3.3) 



(cf. Eq. (25.16) of MTW §). 

Throughout the inspiral and transition regimes, the 
body moves along a nearly circular orbit; its change of 
radius during each circuit around the hole is Ar <C r. 
(Only after the body is well into its final plunge toward 
the hole does Ar become comparable to r.) This near- 
circular motion guarantees that the ratio of the energy 
radiated to angular momentum radiated is equal to the 
body's orbital angular velocity ||: 



dE _ ^dL 
df df 



(3.4) 



Correspondingly, in and near the transition regime, 
which occupies a narrow range of radii around n sco , the 
body's energy and angular momentum are related byn 



E = Eu 



^isco£ 



L = L h 



(3.5) 



By combining Eqs. (|3.5| ) and (3.2), we can regard the 
body's effective potential as a function of f and the dif- 
ference £ = L — L lsco of its orbital angular momentum 
from that of the isco. 



* In reality, finite-mass-ratio effects, incl uding; those dis- 
cussed in the paragraph preceding Eq. (3.11) below, will alter 
these energy-angular-momentum relations by amounts that 
scale as the first and hig her powers of rj. For example, in 
going from Eq. (3.4) to (3.5), there can be an integration 
constant 5E (which scales as rj or some higher power) so 
E = iJi SC o + SE + f2i sco £. In the presence of such effects, 
we redefine n SC o, -Eisco, and Li sco to be the values of these 
parameters at which the ^-corrected V(r,E,L) has a flat in- 
flection point, as in Fig. j], and fii SC o to be the orbital angular 
velocity at this fi sco . Then Eqs. (3.5) remain valid even for 
finite mass ratio r\. 
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FIG. 1. The gradually changing effective potential V(r,£) 
for radial geodesic motion. Each curve is for a particular value 
of £ = L — iisco. As £ decreases due to radiation reaction, 
the body, depicted by the large dot, at first remains at the 
minimum of the effective potential (£1; "adiabatic regime"). 
However, as £ nears zero (at £ ~ £2), the body cannot keep up 
with the rapid inward motion of the minimum; it lags behind 
in a manner described by the transition-regime analysis of 
Sec. III. At £ ~ £5 the effective potential has become so steep 



that radiation reaction is no longer important, the transition 
regime ends, and the body plunges toward the black hole with 
nearly constant energy and angular momentum. 

Figure |l| shows V(f, £) for a sequence of angular mo- 
menta £1, . . . ,£5 around £ = 0. As £ decreases to £ = 0, 
the minimum of the potential flattens out and disappears; 
and just when it is disappearing, the minimum's radius 
r m i n is moving inward at an infinite rate: dr m { n /d^ — > 00 
as £ -► 0. 

In the adiabatic regime of large £, the body sits always 
at the minimum of the effective potential. Its orbit is 
a slowly shrinking circle, guided inward by the motion 
of the minimum. As £ nears zero and the minimum's 
inward speed grows large, the body's inertia prevents it 
from continuing to follow the minimum. The body begins 
to lag behind, as depicted at £ = £ 2 in Fig- 0- This lag 
invalidates the adiabatic inspiral analysis of Sec. || and 
initiates the transition regime. 

As £ continues to decrease, there comes a point (near 
£5 in Fig. |l|) at which the effective potential has be- 
come so steep that its inward force on the body domi- 
nates strongly over radiation reaction. There the transi- 
tion regime ends, and the body begins to plunge inward 
rapidly on a nearly geodesic orbit with nearly constant 
E and L. The objectives of the following subsections 
are to deriv e a set of equations describing the transition 
regime (Sec. Ill B| ) , and show how the transition matches 
smoothly onto the adiabatic regime at large pos itive £ 
and to the plunge regime at large negative £ (Sec. Ill C). 



B. Equation of Motion for Transition Regime 

Throughout the transition regime, because the body 
moves on a nearly circular orbit with radius close to rj sco , 
and because the body's small mass fi = r]M <C M keeps 
its radiation reaction weak, its angular velocity remains 
very close to fti SC o 



^ = ft ~ ft;. 



dt 



(3.6) 



and its proper time ticks at very nearly the standard isco 
rate 



df _ (df \ 
dt \dt J 



1 - 3/r isco + 2a/ffJ c 2 
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1-3/2 



(3.7) 



cf. Eq. (5.4.5a) of 0]. 

This nearly circular motion at f ~ fj SCO produces grav- 
itational waves which carry off angular momentum and 
energy at very nearly the same rate as they would for 
circular geodesic motion at fj s co . T his means that E and 
L evolve in accord with Eqs. (|3.5[), where 



df 



- KTJ 



and 



32-7/3 



l + o/ 



-3/2 



1 - 3/f is 



2a/r\ 



-3/2 
isco 



(3.8) 



(3.9) 



cf. Eqs. (|J), (|J), (§7i]), and Table | It is the smallness 
of 77 = fi/M (e.g., 77 = 1CP 5 for the realistic case of a 
IOMq black hole spiraling into the 10 6 M Q black hole) 
that makes the angular momentum £ evolve very slowly 
and keeps the body in a nearly circular orbit throughout 
the transition regime [cf. the factors of rj that appear in 
Eqs. and Q3.20|)— which with Eqs. ( |3.11| ) and Q3.22| ) 
imply df / df oc 77 J/D .] 

In the transition regime, the body's radial motion is de- 
scribed by the geodesic equation of motion with a radial 
self force per unit massfj] r/F se if inserted on the right-hand 
side: 



df 2 



2 df 



V F sc\l 



(3.10) 



^This radial self force, like the radiation reaction force that 
drives the inspiral, is produced by interaction of the body 
with its own gravitational field — that field having been in- 
fluenced by the black hole's spacetime geometry; see, e.g., 
Ref. Q . The contravariant radial component of the self force, 
with dimensionality restored using r = Mr and r = Mr, 



is (dp r /dr) sc i f = (fid r/dT 2 ) scli 



(fi/M)(d 2 f/df 2 ) s 
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(We write it as r]F se if because its magnitude is propor- 
tional to rj = fx/M.) 

The radial self force r/F sc .\f is nondissipative (since it 
has hardly any radial velocity with which to couple). 
This contrasts with the 0-directed radiation-reaction 
force, which couples to the orbital angular velocity to 
prod uce a shrinkage of the body's angular momentum 
[Eq. fl3.8Q ] and a corresponding decrease of its energy, 
dE/df — Cld^/df. Because the radial force is nondissi- 
pative, it is of little importance. It can be absorbed into 
the nondissipative effective potential term — \dVjdr in 
the equation of motion. Doing so will not change the 
general character of the effective potential, as depicted in 
Fig. ^; it will merely change, by fractional amounts pro- 
portional to 77, the various parameters that characterize 
the effective potential: the location fi SCO of the innermost 
stable circular orbit (at which the £ = effective poten- 
tial curve has its inflection point), the values at the isco 
of the orbital energy and angular momentum E isco and 
Ljsco, and the constant a defined below. There will be 
other 0(77) changes in f; sco , _Ei sco , L- lsco and a caused by 
the body's own perturbation of the hole's spacetime ge- 
ometry Jll|-p^[. In this paper, we shall ignore all such 
changes, and correspondingly we shall neglect the radial 
self force 7?F S eif ■ 

We shall describe the body's location in the transition 
regime by 



R = 



(3.11) 



Throughout the transition regime both R and £ are small, 
and correspondingly the effective potential can be ex- 
panded in powers of R and £. Up through cubic terms in 
R and linear terms in £ (the order needed for our analy- 
sis), the effective potential takes the form 



V(R,0 



2a 



R 3 - 2/?i?£ + constant 



(3.12) 



where a and are positive constants that we shall eval- 
uate below. Note that for £ = 0, this is a simple cubic 
potential with inflection point at R = 0, i.e. at f — fi SCO ; 
and note that for £ > 0, it acquires a maximum and a 
minimum, wh ile fo r £ < it i s mo notonic; cf. Fig. [l| By 
inserting Eq. ( 3.12 ) into Eq. ( 3.1C ), setting f = f isco + R, 



and neglecting the radial self force or absorbing it into 
fi sco , a and (3 as described above, we obtain the following 
radial equation of motion: 



d 2 R 
df 2 



= -aR l 



(3.13) 



By then s ettin g f = at the moment when £ = and 
using Eq. (3.8) for the rate of change of £, so 

(3.14) 



£ 



-TjKT 



we bring our equation of motion into the form 

d2R »2 a - 

dr z 



We shall explore the consequences of this equation of 
motion in the next subsection, but first we shall deduce 
the values of a and (3. 

The constants a and (3 can be evaluated from t he fol - 
lowi ng re lations, which follow directly from ( 3.12 ), ( 3.11 ) 
and (3.5): 



a 



1 / d 3 V(f,E,L) 
4 I 9f3 



1 d 2 V(r,E,L) ^ d 2 V(r,E,L) 

2 \ dldr dEdf 



(3.16) 



(3.17) 



By inserting expression (3.2) into these relations, we ob- 
tain a and (3 in the limit rj = fi/M —> 0: 



ISCO 



(r 2 + 2[a 2 (E 2 - 1) - L 2 ]f + 10(1 - aE) 2 ) 

V / i 
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for a = , 



f3 = -j— ((L — a 2 Ett)f - 3(1 - aE)(l - a£lj) 

f-\ / i 

ISCO 



36V3 



for a = 



(3.18) 



(3.19) 



Here f; sco and S^i sco are given by Eqs. (2.5) and (2.6); and 
Lj sco and i?i SC o are expressed in terms of fi SC o by Eqs. (2.7) 
and ( |2.8| ). Numerical values of a and j3, computed from 
these equations, are tabulated in Table |. 



C. Solution for Motion in the Transition Regime 

Th e equation of motion in the transition regime, Eq. 
( |3.15 ), can be converted into dimensionless form by set- 
ting 



R = v 2/5 R a x 



t T 



where 



R = (/? K ) 2/5 a~ 3/5 , 



(a /3 k) 



~l/5 



(3.20) 
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cf. Table [j]. The resulting dimensionless equation of mo- 
tion is 



d 2 X 
~df 2 



= -X 2 - T 



(3.22) 



We seek the unique solution of this differential equation 
which, at early times T <^ — 1, joins smoothly onto the 
adiabatic inspiral solution of Sec. [n| In that adiabatic 
inspiral, the orbit is the circle at t he m inimum of the 
effective potential of Fig. |l| and Eq. ( [3.12] ), i.e., the circle 
at R = \J /?£/ 'a = \J —{3kt\t j a., which translates into 

(3.15) x = \J — T for adiabatic inspiral near the isco. (3.23) 



1 



1.5 

1 

0.5 
X 
-0.5 
-1 
-1.5 



FIG. 2. Dimensionless orbital radius X as a function of 
dimensionless proper time T for an or bit n ear the isco. Adi- 
abatic Inspiral: The analytic solution (3.23) for adiabatic in- 
spiral outside the isco. Transition: The numerical solution 




to the dimensionless equation of motion (3.22) for the transi- 



tion regi me in the vicinity of the isco. Plunge: The analytic 
solution ( 3.25) for the orbital plunge inside the isco. 




FIG. 3. Same as Fig. ^, but drawn on a different scale. 

We have not been able to find an ana lytic formula for 
the solution to the equation of motion ( 3.22f ) that joins 
smoothly onto this adiabatic solution, but it is easy to 
construct the unique solution numerically. It is plotted in 
Figures || and ||, along with the adiabatic inspiral solution 
( 3.23[ ) and the plunge solution [Eq. ( |3 .25 ) below]. 

The transition solution is well approximated by adia- 
batic inspiral at times T < — 1, but at T > —1 it de- 
viates from adiabatic inspiral and evolves smoothly into 
a plunge. The solution diverges (X — > — oo) at a finite 
time T = T plungc ~ 3.4120 

In the plunge regime, radiation reaction is unimpor- 



*The divergence of X at T = T p iun g c does not imply a 
divergence of r or any other physical quantity. Rather, it 
marks the breakdown of the transition approximation at ve ry 
large values of \X\, \X\ > X brcak ~ rj 
More specifically, when —X > -Xbreak, 



2 / 5 ; [cf. Eq. Q3.20D]. 
higher-order terms 



tant; i.e., the orbit evolves inward with (very nearly) con- 
stant orbital angular momentum L ~ -Zjfinai and energy 
E ~ -Efmai (which we evaluate below); i.e., the orbit is 
well approximated by geode sic fr ee fall. In the dimen- 
sionless equation of motion ( 3.22j ), the free-fall approx- 
imation translates into neglecting the last term, T, so 
d 2 X/dT 2 = —X 2 , which has the analytic first integral 



dX/dT 



constant 



-X 3 . 



(3.24) 



For large \X\, the constant can be neglected and we ob- 
tain the analytic solution 



X 



-G 



(^plunge 2"*) 



for plunge near the isco, (3.25) 



which is plotted in Figs . || and_p. 

Combining Eqs. Q, (flj), and ( |3~20| ), one finds 
that throughout the transition regime, the energy and 
angular-momentum deficits (i.e., the deviations of E and 
L from their isco values) scale as jfl h . In particular, the 
final deficits in the plunge stage are given by 



-^final -^isco — (^^0 -^plunge 
-^filial -^isco ^isco(^7"0-^pl 

where, as was noted above, 

^plunge — 3.412 . 



4/5 



(3.26) 



(3.27) 



IV. GRAVITATIONAL WAVES FROM 
TRANSITION REGIME, AND THEIR 
OBSERVABILITY 

The gravitational waves emitted in the transition 
regime are all near the orbital frequency 27rf2i SCO and its 
harmonics. The strongest waves arc at the second har- 
monic (twice the orbital frequency): 



2tt 



ttM 



(4.1) 



in the Taylor expansion ( 3.12 ) become important and stop 
the divergence. Well before this (in fact, throughout the 
range 1 <C — X <C -Xbreak) both the transition approxima- 
tion and the free-fall approximation (E = -Efl na i = constant, 
L — Lflnai = constant) are valid, so these two approximations 
can be matched in this regime to obtain a solution valid all the 
way down to the horizon. The same type of breakdown also 
occurs at the other asymptotic limit +X > -Xbre ak: The tran- 
sition regime's adiabatic-inspiral equation (3.23) breaks down 
and must be replaced, via matching at 1 <C X <C -Xbreak, by 
the exact Kerr metric's adiabatic inspiral formulae M. 



r. 



We shall compute their properties. 

The transition waves last for a proper time At = 
MAf = M77 _1 / 5 f AT, during which the body spi- 
rals inward through a radial distance Ar = MAR = 
Mrj 2 / 5 R AX, where AT covers the range T ~ —1 to 
~ 2.3 and AX covers the range X ~ 1 to X ~ —5 (Fig. 
1); i-e., 



AT = 3.3 , AX = 6 



(4.2) 



Correspondingly, neglecting any cosmological redshift, 
the duration of the transition waves as seen at Earth 



is 



At = 



M 



(df/dt)i 



-T) 



4 / 5 f„AT 



(4.3) 



and their frequency band is Af = (1 / ttM) (dCt / df)j SC O Af , 
which, using the above expression for Ar and Eq. (2.1) 
for 0(f), gives 



2ttM 



tf sco V^~ori 2/5 RoAX 



(4.4) 



The total number of cycles of these transition waves is 



N r . 



/At 



ir(df/dt)i : 



-1/5 AT 



(4.5) 



These second-harmonic waves arriving at Earth have 

3 (2^ jr fdt + tp+), h x = 



xampCOs(27r J fdt + fx), where ip+ and ip x arc con- 



the form = 
ft 

stant phases. The amplitudes /i+amp and h x amp depend 
on the source's orientation. When one squares and adds 
these amplitudes and then averages over the sky ("(...)"), 
one obtains an rms amplitude 



Here 5Sh(f) is the spectral density of LISA's strain noise 
inverse-averaged over the sky^] and 1/At is the band 
width associated with the waves' duration At. 

The noise spectral density Sh(f) for the current straw- 
man design of LISA has been computed by the LISA 
Mission Definition Team An analytic fit to this 

Sh(f), after averaging over some small-amplitude oscil- 
lations that occur at / > 0.01 Hz, is the following: 



S h (f) = 



(4.6 x 10~ 21 ) 2 + (3.5 x 10^ b ) 



26\2 



+ < 3 - 5 * 10 - 19 > 2 (ts)' 



Hz" 1 . (4.9) 



The rate for fj, ~ 10M Q black holes to spiral into 
M ~ 10 6 M Q black holes in galactic nuclei has been 
estimated by Sigurdsson and Rees their "very con- 
servative" result is ~ one event per year out to 1 Gpc. 
The inspiraling holes arc likely to be in rather eccentric, 
nonequatorial orbits [ fl6| , for which our analysis needs 
to be generalized. If, however, the orbit is circular and 
equatorial and the holes are at 1 Gpc distance, then the 
above formulas give the numbers shown in Table |n|. 

As shown in the table, the signal to noise for this source 
is of order unity. With some luck in the orientation of 
LISA, the orientation of the source, the distance to the 
source, and/or the holes' masses, a S/N of a few might 
occur. Since the signal would already have been detected 
from the much stronger adiabatic inspiral waves, this sig- 
nal strength could be enough to begin to explore the de- 
tails of the transition from inspiral to plunge. 



V. CONCLUSIONS 



.rms 
amp 



{hi 



X amp/ 



a/2 



(4.6) 



which is related to the power being radiated into the 
second harmonic by E 2 = 4ttD 2 (h™ s p ) 2 (2nf) 2 /(32tt); 
cf. Eq. (35.27) of MTW §. Here D is the distance 
to the source. Equating this to the radiated power 



E-2 



(32/5)ry 2 0; 



•<oo,2 



where £oo 2 is a relativis- 



tic correction factor listed on the first line of Table IV of 
[B], we obtain the following expression for the waves' rms 
amplitude 



j rms 
amp 



(4.7) 



The signal to noise ratio S/N that these waves produce 
in LISA depends on the orientations of LISA and the 
source relative to the line of sight between them. When 
one squares S/N and averages over both orientations, 
then takes the square root, one obtains 



N 



u rms 
amp 



y/5S h (f)/At 



(4.8) 



Our analysis of the transition regime has been confined 
to circular, equatorial orbits. This is a serious constraint, 
since there is strong reason to expect that most inspiral- 
ing bodies will be in orbits that are strongly noncircular 
and nonequatorial p6[ . Our estimated signal-to-noise ra- 
tio, S/N ~ 1, for LISA's observations of the transition 



° i.e., l/(5Sfi) = average over the sky of l/(spectral den- 
sity). The factor 5 in this definition is to produce accord 
with the conventional notation for ground-based interferome- 
ters, where Sh{f) denotes the spectral density for waves with 
optimal direction and polarization. In the case of LISA, at 
frequencies above about 0.01 Hz, the beam pattern shows 
sharp frequency-dependent variations with direction due to 
the fact that the interferometer arms are acting as one-pass 
delay lines rather than optical cavities, and this produces a 
more complicated dependence of sensitivity on angle than for 
ground-based interferometers. As a result, (as we have 
defined it) is the spectral density for optimal direction and 
polarization only below about 0.01 Hz, not above. 



i: 



regime from a circular, equatorial orbit at the plausible 
distance ~ 1 Gpc suggests that for more realistic orbits 
the transition regime might be observable. This prospect 
makes it important to generalize our analysis to more 
realistic orbits. 

Full analyses for equatorial, noncircular orbits and for 
nonequatorial, circular orbits can be carried out using 
techniques now in hand: the Teukolsky formalism, and 
computations of the orbital evolution based on the en- 
ergy and angular momentum radiated down the hole and 
off to infinity (see, e.g., Ref. [gj and references therein). 
For nonequatorial, noncircular orbits, the analysis should 
also be possible with existing techniques — up to an un- 
known radiation-reaction-induced rate of evolution of the 
Carter constant. That unknown quantity could be left as 
a parameter in the analysis, to be determined when cur- 
rent research on gravitational radiation reaction ]l7| , p^| ,p[ 
has reached fruition. 

When this paper was in near final form, we became 
aware of a similar analysis, by Buonanno and Damour 
p9| , of the transition from inspiral to plunge. Whereas 
we treat the case of infinitesimal mass ratio rj <ti 1 and 
finite black- hole spin —1 < a < +1, Buananno and 
Damour treat finite 7/ (0 < i] < 1/4) and vanishing spins 
a = 0. Both an alyse s give the same dimensionless equa- 
tion of motion ( 3.22| ) for the transition regime. 
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TABLE I. Dimensionless parameters characterizing the isco and the transition regime of inspiral. The values of £ are from 
numerical solutions of the Teukolsky equation by L. S. Finn (first line of Table II of Ref. B). 



(1 


7"isco 




£ 


a 


P 


K 


Ro 


To 


-0.99 


8.972 


0.03863 


1.240 


0.0001543 


0.006626 


0.005013 


3.129 


45.50 


-0.9 


8.717 


0.04026 


1.233 


0.0001732 


0.007070 


0.005527 


3.117 


43.04 


-0.5 


7.555 


0.04935 


1.197 


0.0003070 


0.009730 


0.008966 


3.048 


32.69 





6.000 


0.06804 


1.143 


0.0007716 


0.01604 


0.01955 


2.925 


21.05 


0.2 


5.329 


0.07998 


1.114 


0.001240 


0.02057 


0.02914 


2.852 


16.82 


0.5 


4.233 


0.1086 


1.053 


0.003115 


0.03270 


0.06291 


2.687 


10.93 


0.8 


2.907 


0.1737 


0.9144 


0.01401 


0.06446 


0.2123 


2.326 


5.539 


0.9 


2.321 


0.2254 


0.7895 


0.03447 


0.09039 


0.4214 


2.041 


3.770 


0.99 


1.454 


0.3644 


0.4148 


0.2234 


0.1289 


1.531 


1.284 


1.867 


0.999 


1.182 


0.4379 


0.2022 


0.5127 


0.09568 


2.594 


0.8551 


1.510 



TABLE II. Properties of the second-harmonic, transition-regime gravitational waves from a /i — WMq black hole spiraling 
into a M = 1O 6 M0 black hole (so r\ — /jl/M = 10~ J ) at r = 1 Gpc distance. The values of £00,2 are from numerical solutions 
of the Teukolsky equation by L. S. Finn (first line of Table IV of Ref. [§). 
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/, Hz 


A/ 
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At, sec 
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£00,2 


L nns 
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(#) . 10 

\ J v / rms 


-0.99 


0.002496 


0.033 


9300 


23 


1.029 


2.0 


1.2 


-0.9 


0.002601 


0.033 


8800 


23 


1.020 


2.0 


1.2 


-0.5 


0.003188 


0.037 


7000 


22 


0.9734 


2.3 


1.4 





0.004396 


0.044 


4800 


21 


0.8957 


2.7 


1.6 


0.2 


0.005167 


0.047 


4100 


21 


0.8535 


2.9 


1.6 


0.5 


0.007016 


0.054 


2900 


21 


0.7653 


3.4 


1.6 


0.8 


0.01123 


0.062 


1900 


22 


0.5914 


4.1 


1.3 


0.9 


0.01457 


0.063 


1700 


24 


0.4617 


4.3 


1.1 


0.99 


0.02354 


0.051 


1800 


43 


0.1656 


3.6 


0.72 


0.999 


0.02829 


0.037 


3400 


96 


0.06128 


2.4 


0.58 
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